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Abstract 



Porous silicon, obtained by electrochemical etching, has been used as a substrate for the growth 
of nanoperforated Nb thin films. The films, deposited by UHV magnetron sputtering on the 
porous Si substrates, inherited their structure made of holes of 5 or 10 nm diameter and of 10 to 
40 nm spacing, which provide an artificial pinning structure. The superconducting properties were 
investigated by transport measurements performed in the presence of magnetic field for different 
film thickness and substrates with different interpore spacing. Perpendicular upper critical fields 
measurements present peculiar features such as a change in the Hc2± {T) curvature and oscillations 
in the field dependence of the superconducting resistive transition width aX, H k. 1 Tesla. This field 
value is much higher than typical matching fields in perforated superconductors, as a consequence 
of the small interpore distance. 



PACS numbers: PACS: 74.25.Fy, 74.78.-w, 74.25.Qt, 81.05.Rm 
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I. INTRODUCTION 



In order to increase the critical parameters of both low- and high-T^ superconductors, in 
the last years great efforts were directed towards the optimization of the vortex confinement 
in these systems. The enhancement of the vortex pinning can be obtained using periodic 



arrays of different types of artificial defects in superconducting thin films, such as periodic 



variation of thin film composition 



l| or heavy ion irradiation |[2|]. Progress in the fabrication 



of nanostructures by lithography has given access to a wide variety of well controlled vortex 
regular pinning arrays at sub- micron length scales 3|, |J, |5|. Experiments performed on 



this kind of samples revealed the presence of peaks in the magnetic field dependence of the 



critical current density Jc{H) 



3, 13, ISl at matching fields Hn = nHi, where n is an integer 



and Hi is the value of the field at which the vortex lattice spacing is equal to the period 



of the regular array of submicrometer defects, a: Hi = 2$o/v^a^. Recently a radically 



different approach to nanostructures fabrication based on self assembled growth 



91], such 



as chemically anisotropic etching of single crystals, attracted much attention also in the 
superconducting field. These processes are useful for generating, on large areas, low-cost 
simple patterns of nanostructures in a single step, assuring a high reproducibility as well. 
Moreover, these techniques give the opportunity to produce regular periodic structures with 
features that can be even an order of magnitude smaller than the ones obtained by time 



consuming e-beam lithography process, leading to templates having features comparable 



with the characteristic superconducting lengths. In particular, the reduced dimensions are 
of great importance in the enhancement of the superconducting pinning properties, since 
smaller interpore spacing implies higher matching fields down to lower temperatures. At 
this purpose transport and magnetic measurements were already performed on Nb thin 
ilms grown on nanoporous AI2O3 templates with pores diameter in the range 25-150 nm 



id . I 111 . I12I . Il3l . Il4l |. Anomalies in the M{H) as well as in the temperature dependence of 
the perpendicular upper critical field Hc2±{T) were observed. 

The aim of this work is to present the superconducting properties of Nb thin films de- 
posited on porous silicon (PS) substrates fabricated by electrochemical etching of monocrys- 
talline Si in an HF solution. Porous silicon is constituted by a network of pores immersed 



in a nanocrystalline matrix ISj]. As a result, PS is characterized by a very large chemically 
reactive internal surface area, which makes this material promising for many technological 
applications in different fields, as for instance micro and optoelectronics, gas sensi ng, and 



biotechnology; recently PS was also applied as a substrate for the nanotubes growth [l5l. Il6|. 
The reason to choose PS as a template for superconducting thin films deposition is twofold. 
On one hand Si based substrates are good candidates for the Nb growth. Moreover, the 
characteristic pore features of these PS substrates are the smallest used in superconducting 
field at this purpose. The porous area of our substrates extends, in fact, on about 1 cm^ and 
it is constituted by pores of mean diameter, 0, tunable between 5 and 10 nm, with interpore 



spacing (the distance between the centers of two nearest pores), A, that can be varied from 
10 to 40 nm. It is worth to under hne that the regularity of the pores arrangement in PS 



substrates is 



ower than the one observed in AI2O3 templates obtained by electrochemical 



oxidation [171 . Il8l . Il9l | . However it has been recently shown that the interaction of the vortex 
lattice with templates of artificial pinning centers presenting a short-range order can also 
give rise to peculiar behavior of both critical magnetic fields and critical currents ^20|] . In the 
same way, thin Nb films deposited on PS substrates can inherit their structure, resulting in 
porous Nb thin films with in plane geometrical dimensions, A and 0, comparable with the 
superconducting coherence length, C,(T). This reduced periodicity could reflect in peculiar 
properties of porous thin films, that could be interesting also in view of possible application. 
For instance tuning the samples geometry (A, 0) we expect to be able to control the super- 
conducting critical temperature, to reach higher matching fields and eventually also higher 
critical currents, down to lower temperatures. The superconducting properties of porous 
Nb thin films have been probed by transport measurements in presence of magnetic fields 
applied in the direction perpendicular to the samples surface. As a consequence of the high 
density of the pore network, the (H,T) phase diagrams present a deviation from the classical 
behavior of Hc2±{T), and the ATc{H) behavior shows fingerprints of the periodic anomaly 
typical of superconducting wire networks. Both these effects appear at fields up to if ~ 1 
Tesla. This value is larger than typical matching fields of periodic pinning arrays obtained 



both by lithographic techniques and by using self-organized AI2O3 templates. 



II. FABRICATION AND CHARACTERIZATION 
A. Substrates 

Porous layers were fabricated by electrochemical anodic etching of p-type, 12 Qcm, and 
n-type, 0.01 Qcm, monocrystalline silicon wafers in 48% water solution of HF at the current 
density of 25 mA/cm^. n-type silicon wafers were anodized with light exposition by 100 W 
tungsten lamp as well without light illumination. The anodization time was chosen in the 
range of 0.5 - 4 min in order to get porous layers with a thickness ranging from 0.5 to 4 
/xm. The integral porosity was estimated by gravimetry to be in the range from 30% up to 
60% ^21|]. For p-type silicon wafers the nominal pore diameter and the interpore spacing are 
respectively = 5 nm and A = 10 nm. In the case of n-type silicon wafers the nominal pore 
diameter is = 10 nm, while the values of A are expected to be 20 nm for light exposed 
substrates and 40 nm for samples obtained without illumination. 

The Electron Back Scatter Diffraction (EBSD) technique has been used to investigate 
the microstructure of both p-type and n-type porous silicon wafers. EBSD patterns were 
obtained from LEO EVO 50 Scanning Electron Microscope with an Oxford Inca detector. 
In order to optimize the contrast in the diffraction pattern and the fraction of electrons 
scattered from the sample, the substrates were tilted so that an angle of 70° is formed 
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between the normal to the sample surface and the electron beam [22]. EBSD patterns were 
acquired on porous and non-porous regions of the two different substrates. Fig. [1] shows 
crystal orientation map of a 3 mm by 0.5 mm area for a n-type substrate with a spatial 
resolution of about 80 /im^. The scanned area was selected crossing the interface between 
the etched (porous) and un-etched (non-porous) silicon wafer. The color uniformity of the 
map reveals that in this area all the crystalline components have (111) orientation. Similarly, 
the same analysis performed on p-type substrate showed (100) crystallographic orientation 
on both porous and non-porous region. In both cases the etching process does not influence 
the crystalline properties of the substrates. 

B. Nb films 

Nb thin films were grown on top of the porous Si substrates in a UHV dc diode magnetron 
sputtering system with a base pressure in the low 10~^ mbar regime and sputtering Argon 
pressure of 1 x 10"'^ mbar. Reference films were grown in each deposition run on standard 
non-porous Si (100) substrates. In order to reduce the possible contamination of the porous 
templates, the substrates were heated at 120°C for one hour in the UHV chamber. The 
deposition was then realized at room temperature after the cool off of the system. The films 
were deposited at typical rates of 0.3 nm/s, controlled by a quartz crystal monitor calibrated 
by low-angle refiectivity measurements. Since the effect of the periodic template would be 



reduced when the film thickness, dm, exceedes the pore diameter, 0, the Nb thickness was 
allowed to range between 8.5 and 15 nm. All the samples deposited on porous Si, their 
names and thicknesses, and the substrate characteristics as well, are summarized in Table 
m For the sake of clarity the samples were named using the initials Si followed by a number 
indicating the nominal substrate interpore distance. A, and by another number for the Nb 
thickness. For example, Si40-Nbl2 is the Nb film 12 nm thick, grown on the porous substrate 
with A = 40 nm. The pore diameter is in this case = 10 nm. 

The morphology of the PS substrates as well as of the Nb thin films deposited on them 
was analyzed by field emission scanning electron microscopy (FESEM). To allow FESEM 
measurements on highly resistive PS substrates, metallization with 8 nm Pd/Au thin film 
was realized. A typical FESEM micrograph (magnification = 500000) of the top of the 
Si20 substrate is shown in Fig. [2|^a). At this high magnification the pores can be clearly 
distinguished as dark areas surrounded by Pd/Au molecular agglomerates. Fig. 2a also 
reveals that the pore entrances are not regular circles. Several micrographs obtained on PS 
templates were analyzed in order to evaluate the pore diameter and the interpore spacing 
distributions resulting from the electrochemical etching. In order to take into account the 
pores shape we evaluated an "equivalent pore diameter" , defined as the diameter of a circle 
having the same surface area of the pore. The latter was obtained by the help of Image Tool, 
a software for image analysis, which enables to supply the surface area of a region drawing a 
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polygon around its outline. The evaluated pore size distribution of the PS layer is reported in 
Figure 2b (mean diameter = 10.35 nm and standard deviation = 1.90 nm). Moreover, 
the mean interpore distance A = 20.42 nm with a standard deviation of a a = 1.80 nm 
(Figure 2c) were estimated. In Fig. [3] a low resolution FESEM image of the Si wafer 
edge is reported. Along this edge, at about 10 fim from the top surface, a change in the 
morphology can be distinguished, indicating the presence of a porous layer of about 10 /im, 
in agreement with sample specifications. All these values are in very good agreement with 
the expected nominal ones and indicate the effectiveness of the porogenic treatment. In Fig. 
Hla secondary electron FESEM image (magnification = 800000) of a dedicated Nb film, 9 nm 
thick, deposited on a Si20 substrate is shown. Again, the pores restrict the solid angle from 
which low energy electrons can escape, thus rendering the pores dark and the top surface 
light. The surface appearance, far from being perfectly ordered, however clearly reveals the 
presence of holes in the Nb layer. 

III. SUPERCONDUCTING PROPERTIES 

The superconducting properties, critical temperatures Tc, and perpendicular upper crit- 
ical fields as a function of the temperature, Hc2±(T), were resistively measured in a ^He 
cryostat using a standard dc four-probe technique on unstructured samples. was taken 
at the 50% of the transition curves. 
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In Fig. [5] are reported the R{T) transition curves of Nb thin films of different thickness 
grown on the porous substrates SilO, namely samples Sil0-Nb8.5 and Sil0-Nbl2, and Si40, 
namely sample Si40-Nb8.5. It can be observed that the critical temperature, T^, decreases 



as the thickness of the Nb film is reduced, as expected [23|. The same Tc{dm) dependence 
is obtained for samples sets deposited on different porous templates. In the same figure the 
resistive transitions of non-porous reference Nb films are also shown. The critical tempera- 
tures of these samples are higher than the of the film with the same thickness deposited on 
the porous substrates. This difference is less evident for the sample with the larger interpore 
distance, since the ratio between the non-porous and the porous area is higher in this case 
so that the critical temperature results less depressed. This Tc reduction could be related, 



for example, to the highly reactive internal porous surface area 15|. From Fig. [5] it also 
emerges that the Tc difference between porous and non-porous samples is reduced as the 
Nb thickness is increased due to the fact that, for larger dm^ the pores are mostly covered 
and their effect on the electronic transport properties is less important. The data show also 
that the R{T) shapes are strongly effected by the pores and the interpore distances dimen- 
sions. In fact, while sample Si40-Nb8.5 shows a rather sharp R{T) curve, sample Sil0-Nb8.5 
presents a broad double step transition. This difference, at these reduced Nb thickness and 
substrate parameters, can be ascribed to the formation of a superconducting wire network 
made of several thin Nb stripes of different width, and hence to the distribution of transition 
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temperatures for the stripes in the entire sample |24j. l25l. l26l. |27|. 

The (H,T) phase diagrams for our samples were obtained performing R{T) measurements 
at fixed magnetic fields and R{H) measurements at fixed temperatures. In general, the 
perpendicular upper critical field of superconducting films of thickness d obeys to a linear 
temperature dependence, 



where ^oll is the parallel Ginzburg-Landau coherence length at T = 0. The temperature 



dependence of ^\\{T) is ^\\{T) = ^on/y 1 — T/Tc. Since the film dimensions in the xy plane 
are larger than ^\\{T), the expression ([T]) is verified in the whole temperature range. 

Before presenting the results obtained for the Nb films deposited on PS, it is worth to 
show the Hc2±{T) curves for two Nb reference films, as reported in Fig. [61 As expected 
the Hc2±{T) behavior is linear for both the samples. A fit of the data with the expression 
([I]), yields a value of the Ginzburg-Landau coherence length at T = 0, ^o|| = 9.5 nm for 
d^b = 8.5 nm and ,^o|| = 8.7 nm for dj^h = 12 nm. Thus the pore diameter dimensions in our 
PS substrates are comparable with the vortex core dimensions at T=0, ~ ,^o||- This means 
;hat each pore can trap only one fluxon, the saturation number being Ug = 0/2^ (T) < 1 
28| . Subsequently multiquanta vortex lattice [Tl] can not be observed in our systems, which 



are then suitable for the study of commensurability effect at low temperatures and high 
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magnetic fields. 

The Hc2±{T) curves obtained for the Nb films deposited on porous Si templates present 
some peculiarities. Fig. [Tl^a) shows the phase diagram, Hc2±{T), of two samples of different 
thickness grown on the porous Si substrate with A = 10 nm and = 5 nm, namely SilO- 
Nb8.5 and Sil0-Nbl2. Similarly, Fig. [7](b) shows Hc2±{T) dependence of two samples of 
different thickness grown on the porous Si substrate with A = 20 nm and = 10 nm, 
namely Si20-Nbl2 and Si20-Nbl5. In both cases it is clearly visible that the introduction of 
the porous array modifies the upper critical field behavior of the Nb thin film. Hc2±{T) is 
characterized, in fact, by a strong non linear temperature dependence near T^, which becomes 
less evident when the Nb thickness exceeds 0. In the sample Sil0-Nbl2, where the condition 
(iATfe = 12 nm > = 10 nm is realized, the upper critical magnetic field is linear in all the 
temperature range, as shown by the linear fit in Fig. [7](a). The comparison between two Nb 
thin films of the same thickness deposited on Si20 and a non-porous substrate is reported in 
the inset of Fig. Efb). Apart from a Tc reduction the porous Nb film is again characterized 
by a non linear Hc2±{T) dependence. However, for these samples the positive curvature of 
Hc2±{T) close to Tc is the only effect we could notice since the extremely small interpore 
spacing of these two analyzed substrates (SilO, Si20) would produce commensurability effects 
only at very large magnetic fields, namely 24 and 6 Tesla, respectively. These fields values 
are not accessible to our experimental set up. On the other hand the matching field should 
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be considerably smaller if Nb films are deposited on a PS substrate with a larger interpore 
distance. For example, for the porous Si template Si40, Hi ^ 1.4 Tesla, a reachable value 
for our experimental set up. For this reason, to study commensurability effects, two films of 
different thickness were deposited on this substrate, namely samples Si40-Nb8.5 and Si40- 
Nbl2. Moreover, in order to investigate a film with an intermediate Nb thickness, the sample 
Si40-Nbl2 was deliberately oxidized obtaining the formation of about 2 nm thick oxide layer. 
This sample, with a resulting suppressed superconducting critical temperature, was named 
Si40-Nbl0. The R{T) transitions measured for different values of the magnetic fields from 
to 2 Tesla are presented in the inset of Fig. [8]^a) for the sample Si40-Nbl0. The resulting 
(H,T) phase diagram is shown in Fig. ^a). Similarly, in Fig. MJo) the Hc2i_{T) dependence 
is shown for the sample Si40-Nbl2. The upper critical fields are again characterized by a 
non linear temperature dependence near Tc. 

Moreover, it is also interesting to note that for these samples the Hc2±{T) dependencies 
present first a positive curvature followed by a change in the concavity, which occurs in both 
cases at ~ 1 Tesla. This effect is more evident if we plot the Hc2± second derivative versus 
temperature as shown in Fig. [8](a),(b), right scale. In both cases dHl2±/dT'^ shows, in fact, 
an abrupt change at if ~ 1 Tesla, where the anomaly in the Hc2±{T) dependence is present. 



Similar features in (H,T) phase diagram were reported for Nb thin films deposited on A 



substrates 
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12| . as well as on triangular array produced by micellar technique 



2O3 



2Q| 
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Another peculiarity that can be ascribed to the reduced dimensionahty of these samples is 
observed plotting the R{T) transition width, AT^ = %{R = 90% Riok)-Tc{R = 10% Rwk), 
as a function of the field. The AT^ dependence, reported in inset of Fig. MJo) for the sample 
Si40-Nbl2 shows, in fact, a minimum again at -ff ~ 1 Tesla. This feature is in general 
expected at the matching field, since dissipative effects are reduced as commensurability 
occurs, resulting in shaper resistive transitions. Both anomalies of the Hc2±{T) and the 
ATc{H) behaviors a.t H ^ 1 Tesla, make us confident to identify this field as the first 
matching field. Hi. Then, it follows that the period of the porous template is a = 48 nm, a 
value in very good agreement with the nominal interpore spacing of this analyzed sample, 
A = 40 nm. The non-monotonic ATc(if) behavior is also reminiscent of the formation 
of a superconducting wire network. In this case the observed oscillations could be a direct 
consequence of fluxoid quantization in a multiply connected superconductor, an extension of 



Little-Parks oscillations of an individual superconducting loops |29| . The expected amplitude 



of the ATc oscillations can be estimated by the expression Tc[^o/aY [30,], which in our case, 
for a = 48 nm, gives ATc = 180 mK, a value close to the one experimentally observed, 
namely AT^ ^120 mK. 

A fit of the data close to Tc with the expression ([T]), yields a value of the Ginzburg-Landau 
coherence length at T = 0, ^oll = 9.4 nm and ^oll = 8.7 nm, respectively. The values of ^oll 



are significantly smaller than the BCS coherence length of Nb, .^o = 39 nm 
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31| , indicating 



that our films are in dirty limit regime with an electron mean free path of / = 1.38 ^qm / 



^0 ~ 3 nm 



321]. The Ginzburg-Landau parameter, n = A(0)/,^o||) can be estimated using 



the expression k = 0.72X^/1 ^ 10, where = 39 nm is the London penetration depth 



of Nb 31|- Ratios of ^ c) \ \/ a ~ 0.2 and A(0)/a ~ 2.0, measured for a = 48 nm, are larger 



than in previous works 
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12 1 on perforated Nb samples, indicating that we are in presence 



of individual vortex pinning 33|. Therefore the porous template can overcome the lattice 



rigidity, being suitable for the vortex lattice pinning in Nb thin films. 



IV. CONCLUSIONS 



Critical temperatures and perpendicular upper critical fields measurements were per- 
formed on Nb thin films sputtered on different PS substrates obtained by electrochemical 
etching of Si in a HF solution. The control of the etching parameter allows to obtain tem- 
plates with pores diameter and interpore distance tunable in the range = 5 — 10 nm and 
A = 10 — 40 nm, respectively. These templates were used as a system of pinning arrays 
in order to obtain the formation of commensurate vortex structures at high matching fields 
down to low temperatures. Anomalies in the Hc2±{T) behavior, as well as in the field de- 
pendence of the R(T) width, were in fact observed at ~ 1 Tesla, which was estimated as 
the first matching field. Hi. These preliminary results are promising, specially considering 
that PS prepared by traditional electrochemical etching exhibits a sponge structure, so that 
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the holes are not a regularly distributed ordered planar pinning structure. Future work 
will focus on critical currents density studies, carried out both by transport measurements 
on patterned samples and by magnetization measurements. Moreover, due to the extreme 
reduced features size of PS templates, these systems look appealing also for the studying of 
two dimensional superconducting wire networks. 

S.L. Prischepa gratefully acknowledges the CNR-Italy for financial support within the 
CNR Short Term Mobility Program. 
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TABLE I: Characteristics of porous Si templates and properties of Nb thin films deposited on such 
substrates. A indicates the interpore spacing, the pore diameter, the Nb thickness and Tc the 
critical temperature of the samples. The interpore spacing and the pore diameter are respectively 
A = 10 nm and = 5 nm for p-type silicon wafers, and A = 20 nm up to 40 nm and = 10 nm 
for n-type silicon wafers. 

Sample (nm) A (nm) dm (nm) Tc (K) 



Sil0-Nb8.5 5±1 10±2 8.5 2.48 

Sil0-Nbl2 5±1 10±2 12 4.68 

Si20-Nbl2 10±2 20±4 12 4.06 

Si20-Nbl5 10zb2 20±4 15 6.18 

Si40-Nb8.5 10±2 40±6 8.5 3.49 

Si40-Nbl0 10±2 40±6 10 5.17 

Si40-Nbl2 10±2 40±6 12 5.94 
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FIG. 1: (color online), (a) Scanning electron microscope image and EBSD orientation map of 
an n-type silicon wafer. The scanned area crosses the interface between the etched (porous) and 
un-etched (non-porous) part of the silicon wafer. The thick solid lines in the blue area indicate 
the boundaries between different regions. The color uniformity in the orientation map (blue area) 
reveals that all of the crystalline components in this region have the same orientation. From the 
color key (b) it emerges that the crystal orientation of the analyzed area has the [111] direction 
parallel to the wafer normal direction, as schematically shown in (c). 
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FIG. 2: (color online), (a) FESEM image (magnification = 500000) of the top of Si20 substrate 
covered with a conductive Pd/Au layer, 8 nm thick. Histogram of the pore diameter (b) and pore 
spacing (c) distributions. 
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FIG. 4: FESEM image of a Nb film 9 nm thick deposited on the PS Si20. 
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FIG. 5: Normalized R{T) transition curves of Nb thin films of different thickness grown on the same 
porous SilO substrate, namely samples Sil0-Nb8.5 and Sil0-Nbl2. For comparison, the resistive 
transitions of non-porous reference Nb films of the same thickness are shown. Also, the R{T) curve 
is presented for the sample Si40-Nb8.5. Riok is the resistance value at T = 10 K. 
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FIG. 6: Perpendicular upper critical field i?c2± vs. temperature of non-porous reference Nb films 
with djvb = 8.5 nm (circles) and d^b = 12 nm (squares). The solid lines represent the linear fit to 
the data. 
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FIG. 7: (a) Superconducting phase diagram (H,T) of Nb thin films of different thickness grown on 
the porous SilO substrate, namely samples Sil0-Nb8.5 and Sil0-Nbl2. The linear fits of the data 
of the samples Sil0-Nb8.5 and Sil0-Nbl2 are shown in figure: near Tc for the sample Sil0-Nb8.5 
and on whole temperature range for the second one. (b) Hc2±iT) dependence of two samples of 
different thickness grown on the porous Si20 template, namely samples Si20-Nbl2 and Si20-Nbl5. 
The linear fits of the data of the samples Si20-Nbl2 and Si20-Nbl5 near Tc are shown in figure. 
The inset shows the comparison among two samples of the same thickness, grown on the porous 
template Si20 and on the non-porous template. 
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FIG. 8: (a) Left scale: Superconducting phase diagram (H,T) of the sample Si40-Nbl0. The solid 
line represents the linear fit to the data near Tc. Right scale: dH'^^JdT^ versus temperature is 
presented. The inset shows resistive transitions of the same sample, measured for different values 
of the magnetic fields, from to 2 Tesla. (b) Left scale: (H,T) phase diagram of the sample 
Si40-Nbl2. The solid line represents the linear fit to the data near Tc. Right scale: dHl^JdT'^ 
versus temperature is shown. The inset presents the R{T) transition width, AT^ = T(.{R = 90% 
Riok)—Tc{R = 10% Riok), as a function of the field. The first matching field is indicated. 
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